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We study the propagation of the scalar modes around a Friedmann-Lemaitre-Robertson- Walker 
universe for general modifications of gravity in the presence of a real scalar field. In general, there will 
be two propagating scalar perturbation fields, which will have in total four degrees of freedom. Two 
of these degrees will have a superluminal propagation — with fc-dependent speed of propagation — 
whereas the other two will have the speed of light. Therefore, the scalar degrees of freedom do not 
modify the general feature of modified gravity models: the appearance of modes whose frequency 
depends on the second power of the modulus of the wave vector. Constraints are given and special 
cases are discussed. 



I. INTRODUCTION 

In [l[ , the scalar cosmological perturbation theory around a general Friedmann-Lemaitre-Robertson- Walker (FLRW) 
universe was studied for those general modified gravity models (MGM) 0, [T|, H, HI, 0, HI whose action is given by 

m 2 r 

s = J d x ^ f(R,Q), (i) 

where f(R, G) is a general function, R is Ricci scalar, and Q is the so-called Gauss-Bonnet term defined by 

G = R 2 - 4R MV R' tv + R^RT"? . (2) 

It was found that, in general, the propagation equation for the metric perturbation was acquiring a term proportional 
to fc 4 , where fc is the length of the comoving wavenumber vector k = \k\. This fc 4 -term implies that, depending on 
the sign of the coefficient in front of the fc 4 -term, either the universe is very unstable on small scales or the short 
wavelength fluctuations propagate with speed larger than the speed of light. The background classical-instabilities, 
due to the presence of a negative fc 4 -term, must be forbidden, whereas there still exists controversy as to if the 
superluminal propagation is compatible or not with experimental observations [^, [ToL [Til [T2I . [l3j . Anyhow, the scalar 
cosmological perturbation for the general MGM exhibits highly non-trivial and interesting features that can be used 
to set constraints on the MGM or may even leave detectable imprints on cosmic structures due to the fc -term. 

These results in [l[ were found for the vacuum only, i.e. no fields other than the gravitational field were included in 
the Lagrangian. An obvious and important extension of this work is to add some matter fields into the MGM action 
and to see how the results we found in the vacuum case, are affected by the existence of matter. This is the main 
purpose of this paper. As a first and simple extension, although useful in early-times cosmology, in this paper, we 
will consider a general MGM together with a single real scalar field tp, whose action is given by 

M 2 f If 
S= l£j ^ X ^~ 9 g ' & ~ 2 / ^tpd^tp. (3) 

If w depends on tp, we can set w to be a constant by a suitable field redefinition. Therefore, we will consider w as a 
constant from the beginning. In this setup, when w = 0, the Lagrangian reduces to the model in vacuum ([T]) 1 . We 
can even remove w, if w > 0, by redefining the scalar field as tp = y/w <p. However, we will leave w as a free parameter 
to enable us to consider those cases for which w < 0. Wc do not assume any particular functional form for f(R, Q, ip). 
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1 When w = 0, <p becomes an auxiliary field. From a variation with respect to tp, we obtain ip as a function of R and Q, i.e. tp = tp(R, Q). 
Removing tp in the action by using this relation yields an action which only contains R and Q. Therefore, obtained action belongs to 

GJ. 
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Although {3} is the basic action we consider, for convenience of the actual analysis, we mainly use a different action, 
equivalent to ©, which is given by 



5 = T? / 

l07T 

Ml 



/(A, a, <p) + (R- X)F(X, a, p) + {Q - <r)£(A, a, p) 



RF(X, a, tp) + £?£(A, a, - 1/(A, a, V ) 



i y d 4 x\/^g to d^tpd^tp, 



where A and <r are auxiliary fields and 

F{X,<r,tp) 



SA' «W) = ^, 

V(X, a, <p) = XF(X, a, v ) + cr£(A, a, p>) - /(A, a, p). 



(4) 



(5) 
(6) 



By the following way, we can verify that the action (|4]) is equivalent to ([3]). By variating S with respect to A and a, 
we have the equations for A and er, which are given by 



(R - \)F\ 
(R - X)F a 



(G - <j)F a = 0, 
{9 - <t)& = 0, 



(7) 
(8) 



where F\ = dF/dX etc. If the combination F\£ a — F% does not vanish, the two equations are independent and A and 
a are given by 



A 

(T 



R, 
Q. 



(9) 
(10) 



Eliminating A and a in the original action by using these results, we find that S in ^ reduces to ([3]) : the equivalence 
of these two actions also holds for the equations of motion. If the combination Fa^ ct — F% vanishes, ([7]) and ([8]) are not 
independent from each other. This does not mean we can not eliminate A and a from the original action. Because ([7|) 
is the condition that the action remains the same under the variation, all the possible pairs of A and a that satisfy 
(J7]) give the same action. We can put any (A, a) we want into S as long as (A, a) are the solutions of ([JJ. Obviously, 
A = R and a = Q are the solution and we can put them into ([4]) to find that it recovers the original action ([3]). 

In the next section, we take the action approach to study the cosmological scalar perturbations. That is, we 
will expand the action ^ to second order in the perturbation and reveal new propagation properties due to the 
inclusion of the scalar field. After the classical analysis of the perturbations, we will also consider quantization of the 
perturbations. 



II. ACTION 



A. Background and Perturbation Variable 



We consider a spacctime which slightly deviates from a FLRW universe. We write the metric as 

ds 2 = -(1 + 2a) At 2 - 2a(t)d t l3&t&x l + a 2 (i)(% + 20% + 28,8^) dx i dx j , (11) 

where a(t) is the scale factor and a,P,<p and 7 represent the metric perturbations. We only consider the scalar type 
perturbation, as both the vector and tensor modes' propagation will not be affected — at linear order — by the presence 
of the scalar field, albeit the modified background dynamical evolution. We also decompose the scalar field as ip + Sip, 
where ip is the background value and 8<p is the perturbation. Since / is a function of R and Q, which must be 
expanded, F and £ defined by ([5]) also acquire perturbations, which we denote as SF and S^. For later convenience, 
we define Z and x by the following equation, 

df 

Z=g^, x = a((3 + a^). (12) 

There are degrees of freedom of choosing S t which is a three dimensional hyper-surface of constant t. Changing 
from Ej to corresponds to a time coordinate transformation: t — > t + T(t,x l ). Under this transformation, the 



perturbation variables transform as 



a = a — 1 , 


(13) 


<p = <p — H 1 , 


(14) 


X = X - 7 1 , 


(15) 


Sip = Sip — ipT, 


(16) 


5~F = SF- FT, 


(17) 


<5£ = ^ - £T, 


(18) 


<fz = - zr. 


(19) 



B. First order action 



To obtain the background equations of motion, we need to expand the action (U) at first order in the perturbation 
variables. We find that, at first order, the action is given by 

M 2 r 

S (1) =Y 2 ~ I d 3 xdta 3 [3(-V + 6HF + 12H 2 F + 24H 4 £ + 24ff 2 #£) </> + 24H(F + 4H 2 £ + 2HH£)j> 



6{F + 4iJ 2 £) 4> - (V + 6HF + \2H 2 F + 72H 2 H£ + 72iJ 4 £) a 



6H(F + AH 2 £)a + ZSip] 



d^xdto- 



' 9 / 

-ip (a + 30) + w(fd(p 



(20) 



The variation with respect to the auxiliary fields S\ and Sa can also contribute to It gives the background 

equations A = R and a = G. In (|20[) . we have already replaced A, a by R, G respectively and eliminated terms 
proportional to <5A and Sa. 

By integrating by parts, we find that the condition = yields the background equations which are given by 



3H 2 = — ( - V - 3HF - 12H 3 £ 
F \2 s 



H = 



4irw 



<P 



2F + 8H£ 



M 2 

wip + 'iwHip — Z 

Ion 



—F + HF — 4H 2 (£ — H£) 
0. 



8ttw 



M 
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(21) 
(22) 
(23) 



C. Second order action 



To obtain the perturbation equations, we need to expand the action up to second order in the perturbation variables. 
We denote the second order action as . Although the calculation is straightforward, the expressions at each step 
of the algebraic calculation become very long. Therefore we shall be content with explaining what we do at each 
calculation step and giving the final expression of the second order action which contains only dynamical fields. 

The field variables we can perturb as free variables are a, \i 4>, &<Pi SF and <5£- Note that <5A and Sa can be written 
as linear combinations of SF, St; and Sip by inverting the relations 

SF = F X S\ + F a Sa + F^Sip, (24) 
St, = F a SX + taSa + Z v 5(p. (25) 

Just for practical convenience, we will use SF and St; instead of SX and Sa. 

By using those fields (a, x, <j), Sip, SF, 8£), we can expand the action up to second order in those fields. Next, 
we chose a suitable gauge to make the calculation easier, the so-called Modified Gravity Models Gauge (MGMG), 
where SF + 4H 2 St; = 0. This gauge condition uniquely fixes the time slicing. By this gauge condition, we can remove 
SF from S^ 2 \ so that is a quadratic functional of {oim, XMi <Am> S<Pmi ^£m}; where the index M indicates that 
we are working in the MGMG. After a few integration by parts, we find that um,Xm and S^m are auxiliary fields. 
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Eliminating these fields by using their equations of motion, becomes a functional of 4>m and S(pM only. Performing 
other integrations by parts, we find that can be finally written as 

M 2 f _ 
S = lfor / dM3a? ^ W { Aab{t ^Vb + a(t)- 2 B(t)e ab VV a ■ W b - a(t)~ 4 D ab (t)AV a AV b - a(t)~ 2 E ab {t)VV a ■ W b 

+ C(t)e ab V a V b + M ab (t)V a V b ]. (26) 
Here V a with (a = 1,2) is defined as 

Vi = 4>m, V 2 = 5<p M . (27) 

The matrices A, D, E and M are given by background quantities. The matrix e ab is the pure antisymmetric matrix 
with element e 12 = 1. Furthermore, the only non-zero matrix element for D is Du. In the short wavelength regime, 
A, B, D and E alone determine the propagation nature of the perturbations. Explicit expressions of these matrices 
are given in the appendix. 

From this action, we can derive some important results regarding the propagation properties of the short wavelength 
modes. First, the determinant of A does not vanish in general. Therefore, there are two independent dynamical 
perturbation fields (hence four degrees of freedom). The reason why we get two dynamical fields is clear. In the 
absence of the scalar field, i.e. in vacuum, it was found that there is only one dynamical field coming from the metric 
perturbation. Because we are now adding a dynamical scalar field which has nothing to do with the gravitational 
sector, it is not surprising that this scalar field brings new degrees of freedom. As an obvious extension, we will get 
2 (TV + 1) degrees of freedom for the scalar perturbations if we add N dynamical and independent scalar fields into 
the action. 

In order not to have ghosts, we require that the two eigenvalues of A ab must be positive. This condition is equivalent 
to the following conditions, 

dctAX), trAX). (28) 

Second important result is that D is not a zero matrix in general. Therefore, the fc 4 -term which has already been 
present in the vacuum case is still present even if the scalar field is added into the theory. Explicit calculations show 
that only Du does not vanish. Therefore, among the two dynamical fields, only one acquires fc -propagation. This, 
we think, is the most important contribution of this paper. In fact, since V\ represents the metric and V 2 describes 
the perturbation of the scalar field, we can say that only the metric perturbation acquires the /c 4 -term behaviour. 
We can understand this result as follows. Since the fc 4 -term is already present in the vacuum case, the appearance of 
the fc 4 -term is solely related to the nature of the modification of gravity. As for the scalar fields, it has a standard 
kinetic term and it is hard to imagine that it would give an additional fc 4 -tcrm contribution from its kinetic term. In 
fact, we find that Du is proportional to the combination F\£, a — F 2 (see appendix), which distinguishes classes of 
modifications of gravity. Therefore, F\£ a — F 2 = is still applicable as a sufficient condition for the absence of the 
k 4 propagation. Below, we will consider a general case where F\£ a — F 2 ^ 0, whereas, in the next section, we will 
consider some special case where, for example, F\^ a — F 2 = 0. 

From (|26p . we can derive the equations of motion for V a . Using a Fourier transformation, the equations of motions 
for the two modes are 

(a 3 A ab V b \ k 2 , , k 2 ■ k 4 k 2 (a 3 Ce ab V b Y 1 

^ + —(a 3 Be ab V b ) + ^Be ab V b + D ab -^V b + E ab — V b + ± 3 " ; + -Ce ab V b - M ab V b - 0. (29) 

In the short wavelength limit (fc large), this vectorial equation can be approximated as 

k 2 ■ k 4 k 2 

A ab V b + -^Be ab V b + D ab ^rV b + E ab ^V b = 0. (30) 
a a 1 a- 

From this equation, we can derive a dispersion relation. By denoting the angular frequency for the mode k as u>, the 
dispersion relation is given by 

/ p lA 1.2 \ 

det l-Lu 2 A + iw^Be + —rD + -^E = 0. (31) 
\ a z a 4 a z J 

This is a quadratic equation in to 2 and it has two solutions. As we discussed before, only two modes among four 
acquire k 4 propagation. From (|31| . to leading order in k, we find that the angular frequency for those modes can be 
written as 

a _ B 2 + tr A tr D - tr(AD) k 4 _ B 2 + A 22 D ll k 4 
U ~ dctA a 4 ~ AnA 22 - A\ 2 a 4 ' ( ' 
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Using the explicit expression for the matrix elements given in the appendix, uj is given by 

3 (F + AHO(Fx + 8H 2 F a + 16H^a) a 4 

It is impressive that this expression, in form, is exactly the same as that in the vacuum case [l[ . However, by using (|24[) . 
we can replace, for example, F\ with different background quantities such as F, F ai F v and ip. After this replacement, 
uj explicitly depends on the scalar fields and is not anymore the same as the vacuum case. There is no clear distinction 
between implicit and explicit dependence of the scalar field on uj 2 . To conclude, the introduction of the scalar field 
does not remove the k 4 propagation but changes the background dependent coefficients in the dispersion relation. 

As discussed in detail in [lj, if the right hand side of (|33[) is negative, the perturbation grows exponentially in time. 
The growth rate is higher for shorter wavelength (proportional to A: 2 ). Therefore, the FLRW universe is very unstable. 
If the RHS is positive, the perturbation propagate as a wave. In this case, the propagation speed (group velocity) is 
linearly dependent on k. Therefore, the propagation speed is in general supcrluminal at small scales. 

As for the dispersion relation for the other two modes, we find that uj, to leading order in k, can written as 

. % . * (34) 

B 2 + A 22 D X1 a 2 a 2 ' 

Therefore, these modes propagate at the speed of light. 



III. SPECIAL CASES 



In the previous section, we have found that the scalar perturbations have four degrees of freedom (two dynamical 
fields). Among them, two modes acquire a k 4 propagation which depends on background quantities, whereas the 
other two propagate at the speed of light. However, there are special cases for which we can not apply this generic 
analysis. This is the case if f(R, Q, <p) satisfies the special relation F\^ a — F 2 = 0, or w = 0, or the universe undergoes 
de-Sitter expansion. In the following, we will consider each case separately. 



A. Case I: F A £ CT -F; = Q,w=£Q 

In this case, both B and D a b vanish identically. Therefore, the dispersion relation (|3T|) simplifies to 

k 2 

2 • 



dct I -wM + —Ej = 0. (35) 

Explicit forms of the matrix elements of A and E are given in the appendix. This equation yields a quartic equation 
for uj, which is given by 

k 2 k 4 

det Auj 4 - [tr A tr £ - tx(AE)} uj 2 + — det E = 0. (36) 

or a 4 

This equation always has real solutions for uj 2 . Explicit forms of det A, det E and tr A tr E — tr(AE) in this case are 
also given in the appendix. 

Eq. (|36[) has two different solutions for uj 2 , in general. Interestingly uj 2 = k 2 /a 2 is not generally a solution of ([56^1 
any longer. This can be seen by substituting this ansatz into ([36]). Under this substitution, the LHS of (|36|) becomes 

k 4 fc 4 trAtr£-tr(AE) k 4 det E _ 4 M 2 P H 2 (F(, a - £F CT ) 2 k 4 
drtA + ^ dct A ~ ~3 nw <p 2 (F + 4H£)(F a +4H 2 £ <J f^ 4 ' 

Assuming the universe is not de Sitter (the de Sitter case will be studied as another special case), this does not vanish 
unless F£cr — £,F a = 0. This happens, for example, when / = f(R, Q) — U (ip). 

If at least one of the solutions for uj 2 is negative, the perturbations are classically unstable. The condition for not 
having ghosts modes (|28|) . and the condition that the two uj 2 are both positive are equivalent to 

trA>0, detA>0, det£>0, tr AtxE - tr(AE) > 0. (38) 
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Just for the purpose of seeing how our findings arc applied to some concrete models, let us consider simple examples 
which have been frequently studied in literature. The first one is [Til . [TEl [lt| 



f(R,g,tp)=f(R,<p). (39) 



In this case, we have £ = £ ff = F a = and 



3jM-wFF 2 

detA= ^ J, trAtr.E-tr(A.E) = 2dctA, det£ = dctA. (40) 

(F + 2HF) 



Therefore, the dispersion relation becomes 



k^ 2 



detA u/- — =0. (41) 



Therefore, assuming det A ^ 0, all the four modes propagate at the speed of light. There are, however, well-known 
exceptions, i.e. GR and the Brans-Dicke theory. In GR, F = 1 and F = and hence det A = 0. This means that one 
of the fields docs not propagate. The perturbation fields in (l26l) are defined on MGMG. In GR, this gauge condition 
is automatically satisfied, so that there still remains some gauge degrees of freedom. By using this remaining gauge, 
we can always set one of the two perturbation fields to be zero. 
For the Brans-Dicke theory, f(R, G,<p) is given by 

f(R,g,^) = ^ 2 R. (42) 

Using a new field defined by ip — <p 2 /8, the action indeed reduces to the original one proposed in [l7j |. 



S 



d 4 x^ HjR - "^d^ipd^ . (43) 



Therefore, for Brans-Dicke theory, F = 2-Ktp 2 /Mp and £ = 0. Then MGMG imposes that V2 = SifM — 0. Therefore, 
only the V\ field is dynamical and An and En determine alone the propagation nature of the perturbations. Applying 
the explicit expressions for An given in the appendix, we find that An is given by 

A u = (2t0 + 3) J * a . (44) 
4(tp + Hp) 

Therefore, a ghost appears when w < —3/2, which agrees with well known results [l8| . Furthermore, it can be shown 
that the propagation speed of the perturbations is the speed of light. 
Another similar case is the Lagrangian with 

f(R,g,<p) = F((p)R + {((p)g-U((p), and w jt 0. (45) 

In this case the number of pro pagating fields is only one, as the MGMG corresponds to setting Sip = 0. This case 
was the only one studied in [lj|, although the authors claimed to have studied f(R, G, p) — fi{<p, R) +$,{p)g — U(ip), 
with general /1, which has, instead, two propagating fields . 
A second example for a special case is given by [20, Hl[ , 

f(R,g,< P ) = R + f 1 (g,<p). (46) 

In this case, we have F = 1, F — and 

l2H 2 i 2 ^w(l + 4Hh 

det A = , (47) 

(1 + 6H0 

tr Atr E - tr(AE) = — * (3H 2 + UH 3 £ + 2H - 2H 2 £ + 8HH£), (48) 

(1 + 6H£) 

det E = ^-w — o (l6H(H 2 + H)£ - AH 2 '( + 3H 2 + 4ff) . (49) 

M p (1 + 6H£) V ' 
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Then the dispersion relation becomes 



LO 2 - 



- 1 



2H 

H 2+ JT 



I6n, 




= 0. 



(50) 



We find that two out of four modes propagate at the speed of light. This result is consistent with the discussion we 
did just after (f37|) . As for the other two modes, the square of the propagation speed is 



= 1 



2H 
H 2 



167TW 



M 2 {1 + 4H£)H 2 



(51) 



In vacuum, we would get 1 + 2H/H 2 which agrees with There, it was found that the universe must be accelerating 
in order to avoid a negative c 2 , furthermore the propagation was becoming superluminal during super-acceleration, 
i.e. H > 0. In the presence of a scalar field, (|5"Tj) implies that, depending on the sign of w (1 + 4ff£), the scalar 
field contribution can shift c 2 s toward cither positive or negative values. Therefore, in principle, it is possible to have 
positive and sub-luminal c 2 even in a decelerating universe. 



B. Case II: non-propagating field — w — 

Let us here discuss the case in which w = 0. This is a special case, which can be understood by noticing that in this 
case ip becomes an auxiliary field. In fact, after integrating it out, one can show that the Lagrangian reduces to the 
vacuum f(R,Q) case already discussed in [l[. However it is also possible to keep the field ip, and expand the action 
at second order in all the fields. In this case, choosing the MGMG gauge, one finds only one propagating field, 
as expected. However, the fc 4 -term is still present, and one finds that for large k, the angular frequency is given by 

2 _ 64 (^F x Z <p + 2F <T F (p ^-Z tp F 2 -Fx^-^F 2 )H 2 fc 4 

UJ — : 7 . ( OZ ) 

3 (F + AH 2 i)[{FxZ v - F 2 ) + 8H 2 {F a Z ip - F^) + l6H^ a Z v - $)] a 4 

This case is, in a sense, discontinuous from the general case, as the degrees of freedom reduce by two, since the field 
5fM becomes itself dependent on the field 4>m- Eq. (|52|) implies that the fc 4 -term will be vanishing, for example, if 
F = F(f) and £ = 



C. Case III: de Sitter— H = 



In this case, since both 5F and 6£ are gauge invariant, we can not impose MGMG in general and (I26|) is not applicable 
anymore. Instead of using the variables involved in the MGMG, we find it convenient to use gauge invariant variables 
$ defined by 

• ~£±i5=* (53, 

and dip. Note that in de Sitter all background quantities except for the scale factor are constants. Using these 
perturbation variables, we find that the second order action can be written as 



S (2 



M 2 
16tt 



d 6 xdt a J 



6F(<I> 2 - a- 2 V<P ■ V$} 



8nw , 2 



M 2 



(Sf 



'■VSip ■ V5ip) 



m\<& 2 - m|<M</J - m\{5ip)' 



where m 2 , m| and m 2 are given by 



m 



2 _ 2F(F - l92H e £ a - 96H 4 F CT - \2H 2 F X ) 
F\ + 8H 2 F a + 16H^a : 
2F(F V + AH 2 ^) 



m 



F x + 8H 2 F a + 16ff 4 £ ff ' 



_ 16H% - F X Z V + F 2 - 16H%Z V + 8H 2 F^ V - 8H 2 Z v F a 



(54) 

(55) 
(56) 
(57) 



3 2(F X + 8H 2 F a + l6H^ a , 

We find that both fields propagate at the speed of light. In order not to have ghosts, the conditions F > and w > 
are required. 
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IV. QUANTIZATION 



Up to here, we have considered only the classical theory of the perturbations. Although we have found non-trivial 
and interesting results at the classical level, nothing prevents us from going to quantize the perturbations. An obvious 
important application of the quantization will be the quantum generation of the curvature perturbations during 
inflation. In this case, inflation may be caused by the modification of gravity or by the scalar field ip (inflaton). 
Then, the observed curvature perturbations will be a mixture of the scalar metric perturbations and the scalar field 
perturbations. Another possible application would be preheating [13, H3| after inflation, where the oscillations of ip (if 
it is an inflaton) induces rapid creation of both V\ and V2 quanta. These applications are interesting by themselves 
and will be discussed elsewhere. In this section, we will provide general procedure to quantize the perturbations. 

Our starting point is the second order action for the perturbations, which is given by Eq. (|26[) . To quantize it, it is 
convenient to write the action in terms of canonical fields. Assuming det A 7^ 0, there always exists a time-dependent 
matrix Z ab {t) such that 



a 3 (t) Z T AZ = I, 



(58) 



where / is the identity matrix. Then the fields defined by ip a = [Z ) ,V b are the canonical fields, as their kinetic 
term becomes 



did 3 a; ^-5 ab ij a ip b 



(59) 



However, as Z depends on the time, we also get time derivative of Z contributions coming from the first two terms 
in (f26|) . This yields additional terms each of which has the same form as one of the last three terms in (|26|) . which 
can be absorbed into a redefinition of the background dependent matrices. Therefore, without loss of generality, we 
can write the second order action like 



S< 2 



dt d d x C 



did 3 a 



-^Sabi'atpb + b(t)e a bVipa • V^b - -d ab (t)A^p a Aip b - -H ab {t)Vip a ■ Vip b 
+ c(t)e ab ipai>b - ^m ab (t)ip a ip b 



(60) 



where the background dependent matrices for each term are related to the matrices in (|26[) by Z. Although we 
can further perform an orthogonal transformation in order to make d ab (t) diagonal — still keeping the kinetic term 
diagonal — we leave it as a general symmetric matrix. 
We can define canonical momenta conjugate to ijj a as 



dC 

7r a = -r-j- = ipa~ be ab Aip b + ce ab ip b . 

We impose the canonical quantization conditions, which are given by 

= [n a {t,x) 7 n b {t,y)} = 0, 



(61) 



(62) 
(63) 



The Heisenberg equations of motion, with the help of the above commutation relations, yield the following evolution 
equations for ip a , 



ijj a - 2be ab Atp b + 2ce ab ip b + d ab A 2 ip b - (H ab - bea^A^ + (m ab + ce ab )ijj b = 0. 



(64) 



According to 24| , it is always possible to write solutions of the equations consistent with the commutation relations 

as 



ip a (t,x) = 



d 3 k 



(27T) 



3/2 



U ab (t, k)a b (k)e tk - s + U: b (t, k)al(k)t 



(65) 
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Here, the time-independent operators a a (k) satisfy the following commutation relations, 

[aa(ki),a b (k 2 )} = [a* (fci), al(k 2 )] = 0, 
[a a (fci),a^(fc 2 )] = SabS(ki - k 2 ). 



(66) 
(67) 



The matrix U ab (t, k) is a collection of c-number time-dependent functions which satisfy both the classical equation of 
motion. 



U ab + 2(bk 2 + c)e ac U cb + d ac k A U cb + (H ac - be ac )k 2 U cb + {m ac + ce ac )U cb = 0, 



and the constraints given by 



UU* T 


- U*U T 


= 0, 


UU* T 


- U*U T 


= -2i{bk' 2 + c)e 


UU* T 


- U*U T 


= il. 



(68) 

(69) 
(70) 
(71) 

These constraints are necessary for the consistency of the commutation relations (f6"6"|) , (|67[) with the canonical com- 
mutation relations ([62]) . (|63|) . By using the evolution equations ([68]) . we can check that these constraints are satisfied 
at any time once they are imposed at one time. A simple example of U and U that satisfy the constraints, is given 
by U = I,U = —\I — {bk 2 + c)e at some initial time. However, there arc infinite other combinations that satisfy the 
same constraints. Which U we should use will depend on the physical situation (initial conditions) we are interested 
in. Here we do not specify the form of U. 

As is well-known, there exists a state |0) which is defined by 



o o (*)|0) = 0. 



(72) 



Then we can construct a particle-basis by applying the creation operators a\(k) to the vacuum state. Once this is 
done, we can calculate expectation values of any operator for the state constructed in this way. For example, the 
"vacuum" two-point function is given by 



(0\Mt,x)Mt,y)\0) 



d 3 k 
(2^7 



U ac (t,km c {t, k)e lk < s ~y\ 



(73) 



Before closing this section, let us finally make a comment. The evolution equations (|64p in Fourier space can be 
seen as the equations of motion for an anisotropic charged harmonic oscillator with time-dependent spring force in a 
two-dimensional plane, in the presence of a time-dependent magnetic field. In fact, ipi and ip2 represent the x and y 
coordinate of the point mass respectively, and the magnetic field whose strength is 2(bk 2 +c) (for unit charge) is parallel 
to the z-axis. Therefore, the quantization of (|60|) is equivalent to the quantum mechanics of the harmonic oscillator 
under a time-dependent magnetic field. For time-independent spring forces in the presence of a time-independent 
magnetic field, [25l . [26| gave the analytic wave function and energy spectrum, which, in terms of the Heisenberg 
picture we consider in this paper, means that analytic solutions of (|68[) were found. For the time-dependent case, 
[27l . [28[ found analytic wave function for an isotropic harmonic oscillator. For the general case, finding the analytic 
solution of |68|) is not an easy task and it is beyond the analysis presented in this paper. 



V. CONCLUSION 



We have studied the scalar type cosmological perturbations for modified gravity models with a single scalar field. 
Followings arc our findings. 

• There are two dynamical perturbation variables, one from the gravity sector and one from the scalar field sector. 
As a result, we have four degrees of freedom. However, there are a few exceptions. In GR, as is well known, 
the scalar mode of gravity does not propagate and we have only one dynamical perturbation variable. Also the 
case w = — no kinetic term for the scalar field — leads to only one propagating degree of freedom, which arises 
because gravity has been modified. 

• Even in the presence of the scalar field, the action, expanded at second order in the fields, contains a term 
proportional to fc 4 . Therefore, inclusion of the matter field (at least in the form of a single scalar field) does 
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not eliminate such term. To be more precise, among the four independent modes, two modes acquire a fc 4 
propagation. The propagation speed of these modes is affected by the presence of the scalar field through the 
modification of the background dynamics due to the scalar field. Depending on the modified gravity model 
and background solution, propagation becomes superluminal or the universe is unstable on small scales. The 
propagation speed of other two modes is equal to the speed of light. 

• If the model satisfies the special condition F\^ a — F% = and w ^ 0, the fc 4 -term disappears from the action. 

• For the special case — F% = 0, the propagation of the short wavelength modes are dominated by the 
/c 2 -term. If the combination F£ a — £F a does not vanish, none of the four modes propagates at the speed of light. 
It is very interesting and highly non-trivial that two modes which propagate at the speed of light in the general 
case no longer propagate at the speed of light in the special case. 

• If the universe is de Sitter, irrespective of the MGM, all the modes propagate at the speed of light. Both F and 
w must be positive to avoid ghosts. 

In Fig. Q] we summarize the classification of the MGM according to the properties of the propagation for the 
perturbation fields. 

2 propagating fields 



F^ a -F z a =0 



no k^term 


k 4 term ***** 


F(tp)R + 
§(<p)G- 
U(<p), 
w^O 


f(R.G),\ 
w=0 I 


1 prop 





F^-F^O 



FIG. 1: Classification of the MGM according to the properties of the propagation for the perturbation fields. The models 
in the inner circle contain only one propagating field. The models in the orange (darker) region have modes which have the 
fc 4 -propagation, whereas those in the yellow (lighter) region do not have such modes. 
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APPENDIX A: MATRIX ELEMENTS A ab , D ab , AND B FOR THE GENERAL CASE 

The matrix elements of A a b are given by 

An = — A ia = — , An-^v,—, (Al) 
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where 

an = (F + 4H£) ( 768 ^ H 8 £ 2 + 384 F a H 6 £ 2 + 384 & H 6 F £ + 64 & <p 2 ^w £ H 5 + 48 F x £ 2 H 4 
V M p 

+ 192 Ffj H 4 F £ + 16 H 4 F£ a <p 2 -Jj w + 48 ^ H 4 ^ 2 + 32 F CT <^ 2 i^u; £ H 3 + 24 F x F £ H 2 

n 167T o 167T - o • o 167T o • o o • o 

+ 8 J? 2 ^ <p 2 —jw + 192 -F CT 2 x7jW <p 2 H 2 H 2 - 192 — T w p 2 £ a F\ H H + 24 F a H F 

9 167T • • o 9 167T \ , , „ 

+ 4F x p 2 j^w£H + 3F x F 2 + FF x p 2 j^wj, (A2) 

ai2 = <p (F + 4 H£) (192 H 6 + 96 £ F CT i/ 4 + 12 £ F A i/ 2 + 32 & FH 5 + 16&F H 4 

Mp 

+ 384 H 3 F a 2 H 2 + l6H 3 F a F -384 H 3 £ a H 2 F x + 8 F a F H 2 + 2F X FH + F F x ), (A3) 
a 22 = 768H 7 Fi£ a + 32H 3 FFF a + 384H 5 F£F a + 32H 4 F 2 F a + 64H 5 FF£ a + 192H 4 F£F a + 1152H 6 £ 2 F a 
+ 384 J ff 6 F££ T - 3072H 5 H 2 £(F x £ a - F 2 ) - 768H 4 H 2 F{F x £ a - F 2 ) + F 2 F X + 64H 6 F 2 £ a + 8H 2 F 2 F a 
+ 4H 2 F 2 F X + AHFFF X + 144iT^ 2 F A + 16H 4 F 2 £ a + 2304# 8 £ 2 £ (T + 48H 3 f£f x + 24H 2 P£f x , (A4) 

where J is defined by 

J = 768H 7 F££ a + 32H 3 FFF a + 384H 5 F£F a + 32H 4 F 2 F a + 64H 5 FF£ a + 192H 4 F£F a + 1152H e £ 2 F a 
+ 38AH e Fi£ a - 3072H 5 H 2 £(F x £ a - F 2 ) - MH 2 H 2 ^wp 2 {F x £ a - F 2 ) - 768H 4 H 2 F{F X ^ - F 2 ) + F 2 F X 

1V1 p 

+ 64H 6 F 2 £ a + 8H 2 F 2 F a + AH 2 F 2 F X + AHFFF X + 144# 4 £ 2 Fa + 16H 4 F 2 £ a + 2304H 8 <f + 48H 3 F£_F X 

+ 2AH 2 P£F X . (A5) 

Because A a b is symmetric matrix, A21 = A\i. 
The quantity B is given by 



B = -MHH 2 -^w£ 2 {F + AH 2 £)^^- T Zi. (A6) 



The matrix element of Dn is given by 



Mp 



Du = -l28aH 2 ^ + 4H ^) F ^~ F "\ (A7) 

and all the other components, as already said, are identically equal to zero. 

Except for E22, Eab appears in the dispersion relation at the leading order only for the special case F x £ a — F 2 = 0. 
Therefore, we give their explicit expressions for the special case. Irrespective of modified gravity models, E22 is given 
by (Ml) . 

E u = — : — [256H 5 £ 3 + 256H 3 H£ 3 ~ 64H 4 ^ 2 | + 48H 4 Fi 2 + 128H 3 F£ 2 + \&H 2 ^wtp 2 i 2 



(F + 2HF + 12H 2 £) 2 V Mp 

MH'HF? + Q4HHF£ Z - 32H Z F^ + 2\H Z FF£ + 8HF£^-wip 2 + 16HF 2 £ 



Mp 



+ 16HFFS, + 3FF 2 - 4F 2 £ + ^wLp 2 F 2 I , (A8) 

M p J 

E12 = —■ —■ f - l&H 3 i£ a ^rWip + 32H 2 HiF tp £ a - 32H 2 HiF a £ v - 4H 2 F&^w<p 

{F + 2HF + 12H 2 £)(F a +AH 2 £ (T )\ Mp y a ^ Mp 

4H£F a ^wip + 8HF&F V - 8HFF a £ ip - FF.^wtp J , (A9) 



A If. ' - Mp 



8nw 



E22 = -TJ2- (A10) 



M p 
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APPENDIX B: MATRIX ELEMENTS OF A ab ,E ab FOR THE SPECIAL CASE 



The matrix elements of A a h are given by 

An= 2(F + AHj) 

(F + 2HF + l2H 2 S,f 



Mp 



(Bl) 



F + 4H£, 16tt 

A V2 = = T—prWip, (B2) 

F + 2HF+12H 2 £M 2 

<B3) 



L P 

As for the matrix elements of E ab -, they are given by (|A8|) . (|A9|) and (jAlOl) . 

We also give expressions for det ^4, det E and tr^4 tv E — tr(AE), which are necessary for deriving the dispersion 
relation. They are given by 

3^w(F + 4H£)(F + AH 2 if 

det A = — ^ , (B4) 

(F + 2HF + 12H 2 £) 

2j$w(F + 4H 2 i) 2 . 

tr AtiE = tr(AE) + £ 2 {l2H 2 FF a + b6H 3 £F a + UH£Fx + SHFF a 

{Fx + 4H 2 F a )(F + 2HF + 12H 2 £,) 

- 8H 2 £F a - 2£F* + 3FF\ + 32HH^F a ), (B5) 



det E = : i r =-^ — UH^F 2 - 2WH\H 2 + H^F'—^wip 



(F + 2HF + 12H 2 ^)" {F x + AH 2 F a Y^ 2 V M p 
128H 3 ZF a Fx^rw<p 2 - 64HH£F x F a ^w<p 2 - 16H£F^w<p 2 - 128H 2 F^F X F C 



p M± p 

o 167T . „ I6n n •■ 167T o o 167T 

M" 2 " W + 32ff ^ F - F ^M2^ - UH 2 FFxF„ W2 - 



QAH^F^w^ 2 - AZ^FFl^wv 2 + 32H 2 ^F a F Xj72 -w^ 2 - 24H 2 FFxF„ J ^wip 
6AH 2 HFF 2 ^w<p 2 + MH 2 F 2 F 2 - ZFF 2 ^wy 2 - 16HFF x F a ^ W <p 2 + 4&*±^v) . (B6) 
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